Radiation reaction at 3.5 post-Newtonian order in effective field theory 
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We derive the radiation reaction forces on a compact binary inspiral through 3.5 order in the 
post-Newtonian expansion using the effective field theory approach. We utilize a recent formulation 
of Hamilton's variational principle that rigorously extends the usual Lagrangian and Hamiltonian 
formalisms to dissipative systems, including the inspiral of a compact binary from the emission of 
gravitational waves. We find agreement with previous results, which thus provides a non-trivial 
confirmation of the extended variational principle. The results from this work nearly complete the 
equations of motion for the generic inspiral of a compact binary with spinning constituents through 
^*», 3.5 post-Newtonian order, as derived entirely with effective field theory, with only the spin-orbit 

corrections to the potential at 3.5 post-Newtonian remaining. 
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The advent of advanced ground-based gravitational wave (GW) interferometer detectors (i.e., advanced LIGO 
and advanced VIRGO) brings an increasing demand for more accurate waveform templates to be used for detecting 
gravitational waves and for extracting information about the parameters associated with a source, such as the masses, 
spins, distance, and sky location. Currently, a goal of the GW source-modeling community is to produce inspiral 
_j. waveforms accurate through at least 3.5 post-Newtonian (PN) order. The PN expansion is a perturbation theory for 
the gravitational field and the binary's motion in the weak- field and slow- motion limits (see [l| for a review). The 
equations for the relative motion of the binary are known already through 3.5PN order, even when including the spin 
angular momenta of the binary's constituents (see, e.g., [J-Q and references therein). However, the gravitational wave 
^ I , flux and the waveform, especially, are not yet known to such a high order for spinning binary inspiral sources. 

High-accuracy waveforms and source-modeling are also important for matching post-Newtonian inspiral waveforms 
| to numerical simulations of binary mergers for purposes of parameter estimation (see e.g., @) and for accurately 
calibrating phenomenological models like Effective One Body Q and hybrid models • These phenomenological 

' models may be used to construct relatively cheap template banks without having to run a prohibitively large number 
. of expensive numerical simulations of binary mergers. This is especially true if combined with an efficient template 
— i ' bank compression and representation scheme such as provided by the Reduced Basis method [l2| . 

The Effective Field Theory (EFT) approach [l3| offers an efficient and algorithmic computational framework com- 
pared to traditional methods [l], [lj] and has rapidly made useful contributions towards the goal of 3.5PN-accurate 
inspiral waveforms. To date, this includes the calculation of the PN corrections to the binding potential, including 
spin angular momenta of the binary's component masses, through 3PN order [2 , Ijjj, [ill. Il5l - [22j . the spinl-spin2 terms 
computed at 4PN [23j |. the leading order radiation reaction force at 2.5PN [24 j, the multipole moments needed for 
calculating the gravitational wave flux through 3PN [25| . and the moments needed for calculating the waveform 
amplitude corrections through 2.5PN [26j . 

In this paper, we calculate the radiation reaction forces that appear at 3.5PN order (spin effects do not enter 
radiation reaction until 4PN). In doing so, we nearly complete the PN equations of motion for the generic inspiral of a 
compact binary with spinning constituents as computed entirely in the EFT framework. Only the spin-orbit correction 
to the potential at 3.5PN order is remaining to be computed in EFT. This is a rather remarkable achievement given 
that the EFT approach was introduced almost eight years ago [l3| . 

Computing radiative effects in the EFT approach presents a unique challenge since the formalism makes heavy use 
of an action formulation of the binary system. More specifically, it is well-known that Lagrangians and Hamiltonians 
are not generally applicable to dissipative systems, which would make computing the PN radiation reaction forces in 
EFT very difficult. In [24j . it was indicated how this might be overcome using a language and notation from quantum 
field theory but was not given a rigorous foundation within a purely classical mechanical context. Nevertheless, as a 
demonstration of the formalism, the 2.5PN radiation reaction force and the gravitational waveform from the binary's 
leading order quadrupole moment were computed [24| and shown to agree with existing results [27l l28j , thus lending 
credibility to the method. Recently, one of us (CRG) gave a rigorous extension of Hamilton's variational principle in 
[29j that yields a Lagrangian (and Hamiltonian) formulation that suitably and correctly describes generally dissipative 
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systems. We use this formalism here, together with EFT, to compute the 3.5PN radiation reaction force and find 
agreement with previously published results [30l . l3lj . This agreement lends non-trivial confirmation for the validity 
of the extended variational principle for dissipative systems described in [2!| (see also the examples given in that 
reference) . 

This paper is organized as follows. Section |TT] gives an overview of the formalisms needed to compute the radiation 
reaction force at 3.5PN order in EFT. Specifically, Section III Al reviews the EFT of compact binary inspirals and 
Section [II Bl reviews the recently formulated extension of Hamilton's variational principle for dissipative systems [29| . 
Section DTI] discusses the computation of the 3.5PN radiation reaction force in EFT where agreement is shown with 
previous results [30l . l3l| . Section IIVI concludes the paper and Append ix [X| outlines in detail the EFT calculation of 
the leading order 2.5PN radiation reaction force of Burke and Thorne 2717 28i| . 



II. OVERVIEW OF EFFECTIVE FIELD THEORY AND DISSIPATIVE MECHANICS 

We begin by reviewing the effective field theory of compact binary inspirals, focusing mainly on the radiation sector, 
and end the section by reviewing how dissipative (e.g., radiative) effects can be handled within the newly developed 
mechanics for dissipative systems (29j . 

A. Effective field theory of compact binary inspirals 

The EFT approach, introduced by Goldberger and Rothstein in separately treats the relevant scales of the 
binary by successively "integrating out" the smaller scales thereby yielding a hierarchy of EFTs that are related to 
each other through so-called matching calculations. The effects of short-distance physics in a large-distance effective 
theory are parameterized in a manner consistent with the symmetries (e.g., general coordinate invariance) as discussed 
in more detail below. 

The slow inspiral of compact binaries has three relevant scales, from smallest to largest: the size of the compact 
objects (COs) R m , their orbital separation r, and the wavelength of the emitted gravitational waves A. The first 
EFT describes the extended masses in the point particle approximation. To incorporate the finite size of the CO 
one appends to the point particle action all possible interaction terms that are consistent with general coordinate 
invariance and reparamctcrizations of the worldline. This leads to a worldline EFT, for one of the COs (e.g., spherical 
and non-spinning), that is described by the action [l3T | 



Sco = -m J dr + C E J dr £ aP 8 ap + C B J dr B al3 B ap + ■■■ (1) 

where £ Q( 3 and B a p are the electric and magnetic parts of the Weyl curvature tensor. The coefficients Ce,Cb, ■ ■ ■ 
are determined via a matching calculation, wherein a chosen quantity is calculated in the effective theory and in the 
long-distance limit of the "full" theory for the CO. One can show that these extra terms in the action contribute to 
the binding potentials (due to induced quadrupolc moments) starting at 5PN for non-spinning COs. We can therefore 
ignore these terms for the work presented here. 

As one "zooms out" from R m to the orbital radius r, the system is described by General Relativity coupled to 
two COs that are each described by the worldline EFT in (TTJ). At this scale, the particles interact with two kinds 
of gravitational perturbations. The first describes the nearly instantaneous potentials that bind the two particles in 
orbit (for further details see [lH). The second is the long- wavelength GWs emitted by the binary. 

As one "zooms out" from r, the binary itself is described in the point particle approximation by a composite object 
[ill, [24| • This composite object radiates gravitational waves from its time-dependent multipole moments that can 
be calculated by matching onto the radiative moments of the binary at the orbital scale. The worldline EFT for this 
radiating object has an action given by [32[ 

5 rad [^,V]= - / drM(r)-^ f dr L ab (r)[n a L b + u^ ab (r)] + ± f dr c«7 a6ai '" a "(r)V ai • ■ • V „£ a6 (O 

^ ^ n=0 ^ 

1 00 r 

+ 2^ j dTc^J aba ^(T)V ai ■ ■ - VaMxn + ■ ■ ■ (2) 



where x m (t) are the worldline coordinates (r beings its proper time), / afca i - a ™ ( r ) anc j jabai—a n ( T ) are t ne symmetric, 
trace- free (STF) mass and current multipole moments, respectively, of the composite object [H, [Hj], f2£ fe is the 
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angular frequency of the body's rotation as measured in a locally flat Lorcntz frame, and uj^ ab are the spin connection 
coefficients, which couple to the total angular momentum L ab = —L ab . Also, lower case roman letters takes values in 
{f , 2, 3}. The mass of the body is taken to be generally time-dependent M(t) as the body may lose rest-mass energy, 
as measured by a distant observer, via gravitational wave emission. A local Lorentz frame is attached to the worldline 
such that Cq(t) = w M (r) = dx^/dr and e^(r) are space-like vectors that rotate with the body and thus account for 
its spin dynamics. The first several Cn' J ^ coefficients are conventionally taken to be 

c o ' = 1' c o ^ = ~~ g' c i ^ = g- (^) 

It is important to note that ([2]) describes any multipolar, extended body that has a size smaller than the wavelength 
of gravitational waves it emits. Therefore, @ is a rather general and model-independent description of such a system. 
However, using the multipole moments computed in the PN expansion for compact binary inspirals [25l [H[ one may 
use @ to study radiative effects, such as radiation reaction, in compact binary systems. 

As with any field theory coupled to point particles, including the EFTs reviewed here, divergences will appear. 

However, there are an infinite number of parameters in the theory (e.g., c£ , cl ) so that divergences can always 
be absorbed into these coupling constants. Interestingly, these parameters can exhibit a classical renormalization 
group flow due to gravitational screening effects, which manifest as logarithmic divergences in the potentials. Unlike 
with traditional approaches for the PN expansion diver gen ces in the EFT approach have a natural place and 
interpretation in the context of renormalization group theory [131 ] . Spin angular momenta for the binary's component 
masses can be included as in Q. 



B. Classical mechanics for dissipative systems 



Determining the evolution of the compact binary in EFT is achieved by integrating out the gravitational perturba- 
tions from the action. In practice, this is accomplished by solving for the gravitational perturbations and substituting 
these solutions into the original action to yield an effective action for the worldlines of the component masses 

mm, 

which is efficiently carried out using Feynman diagrams (l3j . This procedure is well-suited for conservative interac- 
tions, such as for computing PN corrections to the binding potential of a compact binary, but requires modification 
when studying dissipative systems, such as the inspiral of a compact binary due to the emission of gravitational 
radiation, the reason being that dissipative systems generally do not admit Lagrangian or Hamiltonian descriptions. 

Recently, one of us (CRG) introduced a rigorous variational calculus for Hamilton's principle of stationary action 
that includes systems exhibiting dissipative effects and thus generalizes the usual action principle }29j . The derivation 
and the details will not be given here but instead refer the reader to j29j . However, we will review the relevant 
formalism needed specifically for computing the radiation reaction force in this paper. 

The total system formed by the gravitational perturbations h^u and the worldlines of the compact bodies in a 
binary system, xj((t) for K = 1,2, is closed. Only when the gravitational perturbations are integrated out is the 
dynamics of the worldlines open. When integrating out the long wavelength gravitational perturbations at the level of 
the action one must be careful when applying Hamilton's principle of stationary action to the effective action since it 
is formulated by specifying boundary conditions in time, not initial conditions. This observation may seem innocuous 
except that the resulting effective action for the worldlines describes conservative (i.e., non-dissipative) dynamics - 
the Green function for the gravitational perturbations that appear in the effective action are time-symmetric (as these 
are the ones satisfying boundary conditions according to Sturm-Liouville theory) and hence do not account for the 
dissipative effects of radiation reaction [24| . 

To overcome this problem, one formally doubles the degrees of freedom [29[ so that — > (h^i, h^-f) and 
Xk (xki,xk2) then constructs the following action, 

S[x K i,x K2 , h^i, h^ v2 ] = S[x K i, h^ui] - S[x K 2,h^ v2 ] (4) 

where each action on the right side consists of the Einstcin-Hilbcrt action (here, gauge-fixed in the Lorenz gauge and 
expanded around fiat spacetime) and the worldline EFT action $Z§. Integrating out the long wavelength gravitational 
perturbations using Feynman diagrams at the desired PN order (see [2~ij] for the Feynman rules of the radiation EFT 
and how to construct the Feynman diagrams) gives the effective action for the open dynamics of the binary's inspiral. 
After all variations are performed one is free to identify the doubled variables with the physical ones so that, for 
example, xki , xk2 ~> %k ■ This limit will be called the physical limit. 
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It is convenient, but not necessary, to change variables from the (1,2) variables to ± variables where 



xk+ 



x K . 



xki + X K 2 



Xki ~ X K 2 



This change of variable is motivated by the physical limit since then xk- 
physical worldline of the K th body. By computing the following variation [2£ 



and xk- 
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fi[Xl±,X 2 ±\ 



8x K -(t) 



x K ^=0 



(5) 
(6) 

xk where xk is the 



(7) 



one obtains a set of worldline equations of motion that properly incorporates radiation reaction effects [24|, [29| • It is 
important to note that Q receives non-zero contributions from terms in the effective action that are perturbatively 
linear in Xk- or its time derivatives. Therefore, all terms of quadratic order or higher in any "— " variables can 
be dropped from the effective action for the purposes of deriving the worldline equations of motion. We will take 
advantage of this property in the course of the calculations below in Section Mil and Appendix [A] 

In the remaining sections, we compute the terms of the effective action at 3.5PN order by computing the corre- 
sponding Feynman diagrams in the extended action formalism of [2j|. Once the effective action is derived in Section 
III| we then apply the variational principle in <[7|) to obtain the radiation reaction forces on the binary at 3.5PN order. 



III. RADIATION REACTION THROUGH 3.5PN 

In this section we lay out how to calculate the radiation reaction at 3.5PN in the EFT. However, the first non- 
conservative term in the acceleration enters at 2.5PN order and is due to radiation reaction. We include in Appendix 
lAla detailed calculation of the radiation reaction at 2.5PN using the EFT. The 3.5PN terms can be derived similarly 
so we will often refer the reader to Appendix [A] for formulae. 

The relative acceleration a 1 = a\ — a\ for a non-spinning two-body system can be expanded as 



e "2.5 



+ e°a\ + e 



3.5 i 



4 3.5 



(8) 



where e = 1 counts the post-Newtonian order. The motion through 2PN order is conservative, meaning that the 
motion can be characterized by the conserved total energy and angular momentum. Dissipative effects first arise at 
order 2.5PN due to gravitational radiation reaction. The 3PN terms are again conservative while the 3.5PN terms 
are the first post-Newtonian corrections to radiation reaction. 
The first term in ([8|) is just the Newtonian acceleration, 
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«o = (9) 



while the second term is the 1PN correction, derived from the Einstcin-Infcld-Hoffmann (EIH) Lagrangian [3lj 
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where m is the total mass, r\ = m 1 m2/m 2 is the symmetric mass ratio, n l = x 1 jr with 



V) 
(10) 

(11) 

the sep aration 

between the point masses, and v % = dx 1 jdt. The EIH Lagrangian was calculated using the EFT in [l3|, Il5j |. To go 
from Eq. (flT))) to Eq. (TlTj) it is necessary to order-reduce by substituting the C(e°) expression for the acceleration. 
By including the acceleration at C(e 2 ' 5 ) when performing the order-reduction, we get a contribution to the 3.5PN 
acceleration. This will be discussed in more detail below. In the rest of this section we derive the radiation reaction 
through 3.5PN using the EFT. 
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FIG. 1: The Feynman diagram giving the contribution to the radiation reaction from the mass quadrupole. A and B label 
worldlines that have been doubled in the extended variational principle. 



A. Feynman diagrams 

To calculate the radiation reaction, we begin from the worldline action ([2]) and integrate out the electric and 
magnetic parts of the Weyl curvature tensor giving, to the order we require, the effective action 



C-3.5PN _ o , c , c 
^cff — °mq i O mo ~t C q, 



(12) 



where the terms on the right-hand side stand for mass quadrupole (mq), mass octupole (mo), and current quadrupole 
(cq), respectively. As discussed in Appendix [S] at 2.5PN this is accomplished by calculating the diagram in Fig. [T] 
At 2.5PN, each vertex is given by the leading order mass quadrupole moment pi = Iq + 0(e) where 



K 



x K x K 



(13) 



and where K labels the worldlines in the binary K = 1, 2 and the subscript on J iJ denotes the (relative) PN order in 
e. To calculate at 3.5PN, we need to include PN corrections to the mass quadrupole, 



I ij = r o j + elf + 0(e h5 ) 
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m K < x l K x 3 K + e 
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+ 0(e 1 - 5 ) 



TF 



(14) 
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en 



( —(1 - 3?7)z; 2 - -(5 - 8r))— jaV + — (1 - 3?/)(llr 2 i; J ^ - Urrx^)] \ + C(e 15 ), (15) 
\42 7 rJ 21 J J tf 



where TF denotes taking the trace-free part of the enclosed expressions and in the last line we have gone to center-of- 
mass coordinates and fi = rjm is the reduced mass. However, the diagram in Fig. [1] is the same whether one includes 
the 1PN corrections to the mass quadrupole or not, so we get the effective action (|A21[) derived in Appendix [Al 
namely, 



Sm g = jJ dtri(t)igl(t). 



(16) 



where = If — and 1+ = (If + I%)/2, and Ij[ for A = 1,2 is the mass quadrupole moment formed by the 
A th worldline (A labels the doubled worldlines, not the components of the two-body system, which is indexed by 
K = 1,2). By expanding in e, 



-2. 

fJmq ■ — r 



dtlt(t)l^ ]+ (t) 



dtrl(t)4] + (t) 



dtl?_(t)l$ + (t) + 0(e 



1.5n 



(17) 



The first term corresponds to the 2.5PN radiation reaction, calculated in Appendix lAl The second term can be treated 
similar to the 2.5PN term because the same factor of Iq_ will be functionally differentiated when applying Eq. Q. 



In the last term, however, we need to do a bit more work. It is easiest to use Eq. (|14[) and integrate by parts to move 
the derivatives in the last two pieces of the "— " mass quadrupole onto the "+" mass quadrupole in Eq. (|17|) . Once 
that is done, it is straightforward to vary with respect to the "— " coordinates following Eq. (0. This will be done in 
Section HHIJ 

There are two more diagrams that need to be calculated in order to get the effective action necessary to extract the 
3.5PN radiation reaction force. First, there is a contribution from the mass octupole, as shown on the left-hand side 
of Fig. [5J Second, there is a contribution from the current quadrupole shown on the right-hand side of Fig. [5] These 
will be discussed in turn. 
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FIG. 2: The Feynman diagram giving the contributions to the radiation reaction from the mass octupole (left) and the current 
quadrupole (right). 



Evaluating the octupole contribution we have 



iSr, 



1 



2 (6m p ,) 2 



dt I dt'lT(t)(£* k (t)£? m Jt>))l l ^(t>) 



where 



Following the same procedure that is shown in Appendix [51 we can write the integrand as 



/N\ -[Irani. I \ _ rijk , 



^ D (3)AB + 5aD (A)AB + a 2 D (5)AB 



(18) 



(19) 



(20) 



where a is the Synge world function defined in Eq. (|A9|) and the indices are now contracted with Sij . The quantities 
D AB arc given in Eq. (|A7|) and the numerical superscripts on D AB indicate the number of derivatives with respect 
to a. Thus the contribution to the effective action is 



8niG 



dt J dt'lf(t)I ljkB {t') 

dt [ dt' fi k {t)i ijk+ {t') 



>_ D (3)AB + 5aD V)AB + a 2 D (5)AB 

ill + 5*D% + a^DlH 



Changing variables from t' to s = t' — t, holding t fixed, we find using Eq. (|A17I) 



J dtl l l k {t) J ds 



15 dD ret (s) , 15 d 2 D ret {s) 5 d 3 D ret (s) 1 d 5 D ret (s) 



4s 5 ds 



4s 4 ds 2 



4s 3 d3 2 



4s ds 5 



(21) 
(22) 

I ijk+ (t + a). (23) 



Integrating by parts to put the derivatives on Iij k +{t + s) and Laurent expanding to linear order in s, we find 



Smo 



8ttG 



dtrl k {t) / dsD ret {s) 
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_/ iifc+ (i) + _/^ + (t) 



r (4) 



(7) 



(24) 



The first four terms are divergent and are removed by regularization as described in Appendix [A] Keeping only the 
finite term (i.e., the last one) and using Eq. (|A22[) . we obtain 



G_ 
189 



s mo = -^ I dtfl k (t)lj]l + (t). 



(25) 



Since this contribution begins at order 3.5PN we only need the leading order mass octupole moment 7 yfc = I^ k + 0(e) 
with 



T ijk _ 



E 

K 



^k(x^k x ^k x k) stf — _ i^^^) 



STF' 



(26) 



where Sm = mi — mi is the mass difference and STF indicates taking the symmetric, trace-free part of the enclosed 
expression. Varying Eq. (|25|) using Eq. ([7]) gives the mass octupole contribution to the radiation reaction force at 
order 3.5PN, the result of which will be given in Section Hll Bl 

The current quadrupole contribution follows similarly. Evaluating the diagram on the right-hand side of Fig. [2] we 
have 



is = i f-H- 

9 2 \3m p i 



(27) 
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where Bij is the magnetic part of the Weyl tensor, 



Bij = £iaf3\C al3 Jp U X U P 

= e iab R ab jQ + 0(v), 



and the current quadrupole is given by 



K 



(28) 



(29) 



STF 



where L is orbital angular momentum and S is the spin angular momentum. The spin contribution enters at 4PN so 
we can neglect the second term in the current quadrupole above. 
To linear order, the magnetic Weyl tensor is 



Bi 



ab 



{h a O,bj + hbj,a,o) ■ 



Following the same method as discussed in the Appendix, we can write the action as 

2 
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dt / M J%(t)[aD" AB (a)]'j ijB (t') 



dt / dt'ri{t)[oD'; et {a)]'j lJ+ {t') 



Changing variables from t' to s, integrating by parts, and Laurent expanding as before we obtain 



64ttG 



dtTIit) \ dsD ret (s) 



1 7(2) Art , 3 7(4) U\ , 45 7(5) 



Again, the first three terms are divergent and are removed by rcgularization. Therefore, we obtain 

64G 



45 



(30) 

(31) 
(32) 

(33) 

(34) 



B. Radiation reaction force 



We are now ready to vary the action, Eq. (|12[l . to obtain the equations of motion. Since after varying we will set 
all the "— " variables to zero, we just need to vary with respect to the minus coordinates, using Eq. (|7|), or 



dL eS 
dx K i-(t) 



d 
dt 



dL 



eff 



d±Ki-{t) 



x K -=Q 
£k+=£k 



(35) 



where the effective Lagrangian is given by S'off = J dtL c ff(t). 

There are two contributions at order 3.5PN from the mass quadrupole piece, as discussed below Eq. p7|) . The 
easier one is when the lowest order mass quadrupole has the minus coordinates. In that case, we get 



2G 



(36) 



where is the relative coordinate. To get the other mass quadrupole piece from Eq. (|17p . first substitute the 1PN 
part of Eq. (TTJ} into the action for the minus coordinates, integrate by parts to remove the explicit derivatives (so 
that there are no acceleration terms appearing) and then vary with respect to the minus coordinates. This gives 
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The mass octupolc contribution is, from Eq. poll. 



G 5m ■ , ( 7 ) 
63 m J 



while the current quadrupole piece is, from Eq. (|34[) . 

16G 5m 



45 m 



(2xW M jf + x k v^ kl jf - zVe^jf + x 3 x k e m jf 



(38) 



(39) 



There is one more place where we get a contribution to the 3.5PN radiation reaction. We must substitute the 2.5PN 
radiation reaction for the accelerations that appear on the right hand side of the 1PN acceleration in Eq. (|TU|) . Doing 
this gives 



_ 2G 

^reduced r 



(3 + ?])—x k 5 lJ + —i 1 x l x 3 x k + -(1 - iri)v 2 x k 5 t: > + (1 - 3ri)v i x j v k 
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I 
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jk ■ 



Combining all these terms, the 3.5PN acceleration is 
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2^7; fc e lfci jjf } + x fc ^'e tfci jjf } - x 3 v k e k3l J^ + x 3 x k e lkl J^) , 



(5") 

where the all the moments are evaluated at lowest order except for I± A, which is the 1PN contribution. This agrees 
with Eq. (3.15) in Ref. [H|. 



(41) 



IV. CONCLUSION 



In this paper we computed the radiation reaction forces at the 3.5PN order for a compact binary inspiral using 
the effective field theory approach. To accomplish this, we utilized a recently introduced extension of Hamilton's 
variational principle of stationary action [29| that properly incorporates dissipative effects (e.g., radiation reaction) 
within a Lagrangian or Hamiltonian formalism. 

We find agreement between our 3.5PN radiation reaction forces and those of Iyer and Will [30L l3l|. We also derive 
the Burke-Thorne 2.5PN radiation reaction force [13, HH, first demonstrated in the EFT approach in [lH], in detail 
in an appendix. The agreement between our results and previous results gives a strong, non-trivial check that the 
extended action principle formalism for dissipative mechanics (briefly discussed in Section [II Bt its need emphasized 
in [24l ] , and given a proper rigorous framework in [29| ) is the correct formalism for describing dissipative effects in a 
dynamical system. 

Combined with previous results from the EFT community, our work nearly completes the equations of motion for 
a spinning compact binary through 3. 5PN order and derived entirely using the EFT approach. These previous works 
include: the 1PN fl3l [l5j. 2PN and 3PN [13] non-spinning corrections to the Newtonian potential; the 1.5PN 
1 and 2.5PN [H^| spin-orbit corrections; the 2PN @ and 3PN @, i, HH, spin-spin corrections; and the 2.5PN 
24j and now 3.5PN [this paper] radiation reaction forces. The only contribution remaining to be computed through 
3.5PN is from the 3.5PN spin-orbit correction to the potential. 

This work also paves the way for higher order radiation reaction calculations in EFT, including the conservative 
contribution from the radiative back-scattering of gravitational waves off the total mass of the binary spacetime (i.e ., 
the 4PN tail term [3614391 ] ) and the 4PN contribution from the spin-orbit coupling of the current quadrupole [39|, |4(| ■ 
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Appendix A: Radiation reaction at 2.5PN 

To calculate the radiation reaction, we begin by using the worldline action in Eq. ([2]) and "integrate out" the electric 
and magnetic parts of the Weyl curvature tensor. This is accomplished by calculating the diagram in Fig. ((TJ), where 
at order 2.5PN the mass quadrupole moment is given in Eq. ()13l) . In the dissipative mechanics formulation [29j of 
the EFT the contribution to the effective action from this diagram is 



•o2.5PN 
lD cff 



2 (2m pl ) 



dt dt> IJ{t){£ A {t)£ B {t'))l%(t') 



(Al) 



The Feynman rules that tell how to translate Fig. [T]into Eq. (|A1|) are given in 24 1. Here A and B are indices that 
label the doubled variables (see Sec. Ill B|) . which will be called history indices, and the effective action is independent 
of the basis chosen for the fields. For our purposes, it is convenient to work with the basis in which the history indices 
arc A, B = ±. 

The electric part of the Weyl tensor is given in terms of the metric perturbations to linear order by 



£7 r> a /3 

= 77 [didjh 00 (t, x) + dodohijit, x) - d dih j0 (t, x) - d djh M (t, x)} + 0(h 2 ). 



(A2) 



The two-point function in Eq. (|A1[) is then 



(A3) 

where wc have evaluated the radiation fields at the binary's center of mass (taken to be at the origin of our spatial 
coordinates) so that h a p (t) = h a p (t, 0) and have used 



{h i0 (t)hki(t')} cxP i0 « = 0, 



(A4) 



with 



(A5) 



In terms of the field propagators, the two-point function is 

(£#(t)£jg(0) = I {didjdydv + 2P ljkl dld 2 , - i^dldvdv - mdidjdl) D AB (t - t', 0) 
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(Vikdjdodi<do> + rjjididvdk'dw + Tiadjdodk'do* + rjjkdid Q dv d&) D (t - t', 0), (A6) 



where the prime on the spacetime index of a derivative means that it is a derivative with respect to and the matrix 
of (scalar) propagators in the ± basis is 



D AB = 



-iD adv 
-iD ret 



(A7) 



Since the two-point function is contracted with I l J and Ig, which are symmetric and trace- free, a number of these 
terms will drop, leaving 

ri(t)(££(t)£*(t'))I%(t>) = \ri(t) Ud i d j d«d v +2 nik r ljl dldZ)D AB {t-t'M /«(*') 



(A8) 



The propagators can be written in terms of Synge's world function 



a(x<*,x"*) = \ VaP {x a - x' a ){ X P - x'P) = \{t - t'f \{x-x'f 



(A9) 



10 



Therefore, by use of the chain rule, we can simplify Eq. (|A8[) to 



-D" AB + 2(t - t') 2 D^ AB + -(t- t') 4 D^ AB 
2 4 



-D 



/MB 



2 n»AB\ 



(A10) 
(All) 



where the derivatives on D are with respect to er. Plugging this back into the effective action, Eq. (|A1|) . we have 



;r2.5PN 
'off 



dt / dt' rj>(t)i ijB (t') 



(A12) 



Summing over the history indices A, B — ± and using Eq. (| A7|) gives 



2.5PN 
eff 



8m 2 



<# / dt 1 \ri{t)i l0+ (t') 



-l D '; et + (a 2 D'; et )" 



_ n" 

— 2 ad 



2 n " 



7riG / (it / dt'?i(t)I ij+ (t') 



(A13) 
(A14) 



where in the last line we have used m pl 2 = 32ttG and the identity between the retarded and advanced propagators 

Dadv{x, x') = D ret (x', x). (A15) 

The next step is to write the integral in the effective action in a quasi- local expansion about the point t' = i, which 
is when the retarded propagator has non-vanishing support. To do so, define the time difference as 



t' - t, 



(A16) 



so that 



*=-. (A17) 

By changing variables to t and s and by integrating by parts to remove the derivatives on the propagators, we can 
write the effective action as 



S c 2 ff 5PN = SirG I dtl l l(t) \ dsD ret ( S ) 



The retarded propagator is proportional to 



3^ f l ij+ (t + s) \ 3d^_ f l ij+ (t + s) \ ld?_ f l l0+ {t + s) 
' Ads 2 



Ads 
1 

F 4 



2ds 3 



1 d% j+ (t + s) 



ds 4 



D ret (s,6) oc 

s 

so we Laurent expand the terms in the square brackets up O(s), 



-<2.5PN 



8nG J dtl l l{t) J ds D ret (s) 



1 



r(5) 



+ W 32 



10 w 



(A18) 



(A19) 



(A20) 



The hrst three terms are power divergent and are thus zero in dimensional regularization. If we picked a different 
regularization scheme, these terms will end up renormalizing either the mass or other (possibly gauge-violating) terms 
in the action (see [4lj | for further discussion on this point). We thus can remove them and keep only the finite term 
(the last one), 



rr2.5PN 
'-'off 



AttG 



G 



dtlH(t)I§l(t) I dssD ret (s) 



dtlH(t)lgl(t), 



(A21) 
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where we have used 



J dssD ret {s) = i-. (A22) 



To get the equations of motion once we have the above action, we vary in the usual way with respect to the worldlinc 
coordinates for each body. At the end of the calculation we set the history indices to be the same so that all the 
"— " variables are set to zero. Therefore, to get a non-zero result, we just need to vary with respect to the minus 
coordinate as in Eq. ([7]), or 



dL 2 J PN d ( 8L 2 J PN 



dx K i-(t) dt \dx K i-(t) 



_ n (A23) 

£k+=£k 



where the effective Lagrangian is given by S^ 5PN = J dtL 2 ^ PN (t). Using the leading PN mass quadrupole moment, 
Eq. CE), we get 

_2Gm K d?I"(t) ( . 

— 5 dt ' (A24) 

We have been using the mostly-minus convention for the metric, so for the spatial components if J = —8 13 . Switching 
now to contracting indices with 5 %3 we have 

We thus recover the Burkc-Thorne force jHUll using the EFT [H[. 
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